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ABSTRACT 


This  report  considers  the  numerical  determination  of 
the  current  induced  on  a  bent  rectangular  plate  illuminated 
by  an  incident  plane  wave.  Drawing  heavily  on  techniques 
which  are  first  proven  effective  for  the  two-dimensional 
problems  of  scattering  of  TE  and  TM  plane  waves  incident  on 
infinitely  long  bent  strips,  the  authors  develop  a  unique 
staggered  subdomain  scheme  for  representing  the  currents  on 
a  bent  plate.  The  resulting  numerical  techniques  are  simple, 
efficient,  and  highly  convergent.  Computed  current  dis¬ 
tributions  are  given  for  two  frequencies  and  for  several 
angles  of  incidence  for  both  a  flat  square  plate  and  a  bent 
square  plate. 
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INTRODUCTION 


This  report  considers  the  numerical  solution  of  scat¬ 
tering  by  a  bent  rectangular  plate  configuration-  The 
formulation  is  based  on  the  electric  field  integral  equation 
for  the  induced  current  on  the  plate  structure.  The  method 
of  moments  [1)  is  used  to  obtain  a  numerical  solution  for 
the  current,  from  which  all  other  scattered  field  quantities 
can  be  derived  - 

Previous  approaches  to  solving  for  currents  on  flat 
rectangular  plate  structures  [2]  or  their  Babinet  equiva¬ 
lents,  rectangular  apertures  in  an  infinite  ground  plane 
[31,  suffered  from  slow  convergence  and  from  extremely 
complex  formulations.  During  the  course  of  this  study,  it 
was  found  that  most  of  these  difficulties  can  be  overcome 
by  a  judicious  choice  of  both  the  current  representation  and 
the  numerical  treatment  of  the  derivatives  in  the  electric 
field  equation.  The  resulting  technique  is  simple  and 
efficient  to  apply  and  appears  to  be  generaiizable ,  with 
these  same  advantages,  to  treating  arbitrary  surfaces. 

It  may  at  first  seem  unusual  in  a  report  of  this  nature 
to  include  two  chapters  on  scattering  by  infinite  strips  with 
both  TE  and  TM  illumination.  However,  during  the  course  of 
this  investigation,  the  authors  found  themselves  repeatedly 


returning  to  these  problems  to  test  their  hypotheses  and  to 
compare  various  numerical  treatments.  Basically,  the  two 
strip  problems  contain  all  the  same  difficulties  as  the  plate 
problem  except  for  the  fact  that  the  two  problems  are  un¬ 
coupled.  Thus  the  numerical  treatment  of  c he  plate  edges, 
the  bend,  and  the  handling  of  the  derivatives  in  the  electric 
field  integral  equation  have  all  resulted  from  extensive 
testing  on  strips,  where  the  flexibility  in  the  number  of 
unknowns  that  can  be  used  allows  one  to  easily  check  con¬ 
vergence,  dependence  on  subdomain  size,  and  so  on.  Space 
prevents  us  from  documenting  many  of  these  "test  runs,"  but 
it  can  be  said  that  the  sample  results  shown  in  Sections  II 
and  III  are  quite  representative  of  results  obtained  for 
many  different  combinations  of  frequency,  bend  location, 
bend  angle,  and  angle  of  incidence. 

Such  extensive  numerical  experimentation  for  the  plate 
has  not  been  possible  because  of  computer  storage  and  time 
limitations.  However,  when  it  has  been  possible  to  do 
limited  testing,  the  proposed  numerical  method  has  been 
observed  to  be  extremely  stable,  quickly  convergent,  and  tc 
converge  to  currents  which  are  physically  reasonable. 

Section  IV  discusses  the  numerical  formulation  for  the 
bent  plate  structure  and  Section  V  shows  extensive  results 
for  the  current  distributions  on  a  flat  and  bent  plate  for 
frequencies  and  several  different  angles  of  incidence. 
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II. 


NUMERICAL  SOLUTION  PROCEDURE  FOR  A  BENT  STRIP 
KITH  TRANSVERSE  MAGNETIC  ILLUMINATION 


In  this  section  an  integral  equation  is  formulated  and 
the  moment  method  is  applied  to  obtain  a  numerical  solution 
for  the  problem  of  TM  plane  wave  scattering  by  a  bent  strip. 
Effects  of  the  edge  conditions,  the  bend,  and  abrupt  changes 
in  subdomain  size  on  the  numerical  results  are  discussed.  Some 
representative  data  for  currents  on  strips  of  various  config¬ 
urations  and  illuminations  are  given. 

Formulation  of  the  Integral  Equation 

Consider  a  conducting  strip,  infinite  in  the  z-dimension 
and  bent  at  an  angle  a,  as  shown  in  Figure  2.1.  A  plane  wave 
with  only  a  z-component  of  electric  field  is  incident  in  tie 
x-y  plane.  The  incident  field  is  thus  transverse  magnetic 
(TM)  to  the  z-axis.  The  corresponding  scattered  field  is  also 
TM  and  the  surface  current  induced  on  the  strip  flows  only  in 
the  z-direction.  The  surface  current  produces  a  scattered 
vector  potential 

A  (p)  =  ~  |J  (s')  H<2) (kjp-o’ j)ds’  (2.1) 

I  2  0 

strip 

where  u  is  the  permeability  of  free  space,  k  is  the  wave- 

number,  is  the  z-directed  surface  current  density  on  the 

(2  ) 

strip,  and  is  the  Hankel  function  of  second  kind,  zero 


order.  Requiring  the  z-component  of  the  electric  field  to 
vanish  on  the  strip  surface  results  in  the  integral  equation 


_mc ,-. 

■E  ^o)  = 
z 


-juAz(p)  ,  P  =  ss 


(2.2) 


For  a  numerical  solution  of  (2.2),  we  define  unit  pulse 
functions  along  the  portion  of  the  strip  on  the  x-axis  and 
unit  pulse  functions  along  the  c-axis  portion.  These  sets  of 
pulse  functions  are  of  widths 


Ax  = 


(2.3a) 


4c  =  if  ’ 

c 


(2.3b) 


respectively.  There  is  a  total  of  N =  N  +N  pulse  functions. 

x  c 

Thus  the  current  may  be  approximately  represented  in  the  puls 
basis  as  (Figure  2.2) 


N 

X 


J  (p)  =  7  j  p  (p) 

z  I  n  n 


(2.4) 


where 


Pn(P) 


£  s  1  s  + 


(2.3) 


0  ,  otherwise 


and 


r 


s  <  0 
m 


^  =  < 
sc  ' 


xc 


•«l>  ,  S  >  0 

Tcc  m 


(2.8b) 


The  functions  on  the  right  hand  side  of  (2.8)  are  defined  as 


s 

/ 


(2) 


’^pq(ss;  si»  S2^  ~  I  H0  ^kRpq)ds’  ’  ^2*9) 


wh  ere 


p  =  x  or  c . 
q  =  x  or  c , 


R  =  j 

XX 

x-x '  | 

(2.10a) 

R 

xc 

(x  -  C* 

cosa) 

2  5  ■> 

+  c,4sin~a 

r 

(2.10b) 

R 

cx 

[cc  cosa 

-  X  1 

2  2  2  "1 
)  +  c  sin  al 

! 15 

. 

(2.10c) 

R  =  j 

c-c*  ! 

(2 . lOd) 

cc 


Equations  (2.6)  and  (2.7)  constitute  a  linear  system  of  equa¬ 
tions  which  can  be  succinctly  written  in  matrix  form  as 


ZI  =  V 


(2.11) 


where 
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■<£**  •^r»'~-  CF?^ fcsi^A**  V  V-  1  J£p^'.  i 
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V  = 


_mc  , ~  ..  . 

■E  (ss,  )Ax 
z  1 


„inc  ~  .  . 

■Ez  (ssN  )Ax 

x 


an:.*  ., 

Ez  (sSN  +l)Ac 

•  X 


•Einc (ss  ) Ac 
z  N 


(2.12) 


£! 


I  = 


N 


(2.13) 


and  the  elements  of  Z  can  be  determined  from  (2.6)  and  (2.7). 
Solution  of  the  simultaneous  equations  (2.11)  determines  the 
coefficients  of  the  basis  representation  (2.4). 


Effect  of  the  Edge  Condition  on  Numerical  Solutions 

The  behavior  of  current  flow  near  a  conducting  edge  is 
determined  by  the  so-called  edge  condition  [4  ].  Edge  con¬ 
ditions  for  a  bent  strip  and  plate  are  summarized  in  Appendix 
A.  In  the  TM  case,  the  current  flow  is  parallel  to  the  strip 
edges  and  the  edge  condition  requires  that  the  current  be 
singular  there.  Near  the  left  and  right  hand  edges  of  the 


—  U 

strip  in  Figure  2 . 1,  the  current  varies  asjL^  +  si  and  |Lc~sj  ", 
respectively.  At  lower  frequencies,  this  edge  behavior  is 
dominant  over  the  entire  strip,  while  at  higher  frequencies, 
the  behavior  is  confined  more  closely  to  the  edges  of  the  strip 


-U 
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as  Figures  2.3  and  2.4  illustrate  f«r  a  flat  strip  {u  -  G).  I«« 

[  3  ],  it  is  shown  that  the  pulse  basis  set  of  Figure  2.2  enjoys 
a  number  of  advantages  in  a  numerical  determination  of  the  cur¬ 
rent-  First,  the  match  points  are  logically  chosen  at  the 
center  of  the  pulses  so  that  fields  are  not  evaluated  directly 
at  the  edges.  In  other  schemes  vher^  match  points  are  located 
at  the  edges  (such  as  with  half-pulses  or  half-triangles  at 
the  edges),  it  has  been  found  that  anomalous  behavior  in  the 
current  distribution  at  the  edges  results  and  that  the  anom¬ 
alies  in  such  cases  can  be  alleviated  only  by  permitting  the 
edge  current  expansion  functions  to  contain  the  correct  edge 
singularity  [  5  ].  it  seems  to  be  a  general  rule  with  pulse  basis 
functions  that  one  should  always  avoid  computing  the  vector 
(scalar)  potential  at  a  point  where  the  associated  current 
(charge)  is  known  to  se  singular.  In  a  pulse  basis  solution. 


the  pulses  at  the  edge  tend  to  represent  numerically  the  cor¬ 
rect  average  current  contained  in  the  edge  subdomains.  While 


I 


m 

is 


the  magnitude  of  the  average  current  in  edge  subdomains  is 
typically  larger  than  the  actual  current  at  the  midpoint  of 
the  edge  pulses,  the  numerically  derived  current  at  the  edge 
can  be  corrected  a  posteriori  (  5]. 


m 


Effect  of  the  Bend  on  Numerical  Solutions 


The  bend  in  the  strip  also  causes  singularities  in  the 
TM  surface  currents.  As  discussed  in  Appendix  B,  the  surface 


current  found  in  a  numerical  solution  is  the  sum  of  the  actual 
surface  currents  which  flow  on  opposite  sides  of  the  strip. 
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MRONITUOE;  H  -  REAL;  O  -  IMAGINARY 


Induced  current  on  n  ] . OX  strip  illuminated 
Incident  TM  plane  wave. 


In  the  neighborhood  of  the  bend  the  surface  current  must  be¬ 
have  like  that  on  an  infinite  wedge  of  the  same  wedge  angle; 


but  this  behavior  is  different  for  an  interior  and  an  exterior 
wedge  angle.  Only  the  exterior  wedge  angle  causes  singular 
currents,  however,  so  that  the  singularity  in  the  current  can 
be  determined  from  the  exterior  wedge  angle.  In  terms  of 
Figure  2.1,  the  current  in  the  neighborhood  of  the  edge  be¬ 
haves  like  (see  Appendix  A) 


J“L 


J  (s) 


s+0 


Ais  | 


■i+  a 


(2.14) 


While  a  detailed  investigation  to  check  that  (2.14)  is 
satisfied  in  a  numerical  sense  has  not  teen  carried  out,  it 
is  verified  in  Figures  2.3,  2.6,  and  2.7  that  increasing  the  bend 
angle  does  indeed  cause  the  singularity  to  become  more  pro¬ 
nounced  as  (2.14)  suggests.  Furthermore,  it  is  seen  that  the 
behavior  on  either  side  of  the  bend  can  be  very  different. 

For  this  reason  and  for  reasons  similar  to  those  mentioned  in 
connection  with  edges,  it  appears  that  the  pulses  should  indeed 
be  placed  on  either  side  of  the  bend  with  match  points  at 
their  centers,  away  from  the  bend  itself.  Splitting  a  pulse 
around  the  bend  would  be  expected  to  cause  the  same  anomalies 
in  the  current  distribution  as  is  found  when  match  points  are 
placed  at  edges. 


Effect  of  Differing  Subdomain  Sizes  on  Numerical  Solutions 
In  many  numeri cal  procedures,  it  has  been  observed  that 
anomalies  in  the  current  distribution  result  when  abrupt  changes 
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by  a  TM  plane  wave  normally  incident  to  the  section  of  the  strip 


lNaaano 


in  subdomain  sizes  are  made.  In  strip  problems,  enough  subdo¬ 
mains  can  be  used  such  that  if  unequal  subdomain  sizes  are 
necessary  on  opposite  sides  of  the  bend,  the  differences  in 
subdomain  sizes  can  be  made  rather  smhll.  In  the  plate  problem 
to  be  considered,  however,  where  the  number  of  subdomains  in 
each  dimension  is  severely  limited,  it  may  turn  out  that  in 
order  to  adequately  sample  currents  on  opposite  sides  of  the 
bend,  drastic  differences  in  subdomain  sizes  must  be  used.  To 
investigate  the  dependence  of  the  numerical  solution  on  abrupt 
changes  in  subdomain  sizes.  Figures  2.8  and  2.9  show  the  effect 
on  the  solution  of  changes  in  subdomain  size  on  a  flat  strip 
where  the  change  is  made  in  the  central  portion  and  near  an 
edge  of  the  strip,  respectively.  It  is  seen  that  only  a  very 
small  anomaly  occurs  when  the  change  is  made  near  an  edge, 
where  the  slope  of  the  true  current  is  the  greatest.  This 
anomaly  is  so  small  as  to  be  of  no  concern  in  a  numerical 
solution.  This  independence  of  subdomain  size  has  been  veri¬ 
fied  for  a  wide  range  of  strip  sizes  and  incidence  angles, 
provided  that  one  always  begins  with  a  sufficient  number  of 
pulses  (and  match  points)  that  the  variation  of  the  incident 
field  is  adequately  sampled.  Figure  2.10  shows  that  changing 
subdomain  size  at  the  bend  also  has  no  appreciable  effect, 
except  that  the  currents  next  to  the  bend,  which  approximate 
the  average  currant  in  their  respective  subdomains,  always 
appear  to  be  slightly  too  large  in  magnitude. 
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III.  NUMERICAL  SOLUTION  PROCEDURE  FOR  A  BENT  STRIP 
WITH  TRANSVERSE  ELECTRIC  ILLUMINATION 

In  this  section  an  integral  equation  is  formulated 
and  the  moment  method  is  applied  to  obtain  a  numerical 
solution  for  the  problem  of  TE  plane  wave  scattering  by  a 
bent  strip.  The  placement  of  the  current  subdomains  at  the 
edges  and  the  effects  of  the  bend  and  abrupt  changes  in  sub¬ 
domain  size  on  the  numerical  results  are  discussed.  Some 
representative  data  for  TE  currents  on  strips  of  various 
configurations  and  illuminations  are  given. 

Formulation  of  the  Integral  Equation 

Consider  a  conducting  strip,  infinite  in  the  z-dimension 
and  bent  at  ar.  angle  o,  as  shown  in  Figure  3.1.  A  plane 
wave  with  only  a  z-component  of  magnetic  field  is  incident 
along  the  x-y  plane.  The  incident  field  is  thus  transverse 
electric  (TE)  tc  the  z-axis.  The  corresponding  scattered 
field  is  also  TE  ana  the  surface  current  induced  on  the  strip 
flows  only  in  the  s-direction.  Since  the  total  tangential 
electric  field  must  vanish  on  the  strip  surface,  we  have 

-E  ( 0 )  •  s  =  I  -  j u;A ( p )  -  ?$(p)  *s  ,  (3.1) 

5  =  ss 
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Figure  3.1.  Strip  configuration  for  TE  illumination 
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The  vector  and  scalar  potentials  in  (3*1)  are  defined  by 


acB)  =  Jjj*j(p’)  h^2)  (klB-5’ |)ds’ 


(3.2) 


strip 


and 


4  i  e 


f 


0(p)  =  TTT  jp(p’)  Hq2^ (k | p-p *  j  )  ds  ’ 


(3.3) 


strip 


respectively,  where  J(p)  =  s  J(p)  is  the  s-directed  current 
density  and  P  is  the  surface  charge  density  which  is  related 
to  the  current  through  the  equation  of  continuity. 


p‘"p>  -  i  ^sr1 


(3.4) 


To  effect  a  numerical  solution  of  (3.1),  we  divide  both 
straight  segments  of  the  strip  into  subdomains  as  in 
Figure  3.2.  Pulse  functions  representing  the  current  are 
defined  such  that  half  pulses  represent  the  current  on 
either  side  of  the  bend  and  at  the  edges  where  the  current 
vanishes.  The  half  pulses  at  the  edges  improve  the  conver¬ 
gence  of  the  solution  (see  the  following  subsection)  while 
the  half  pulses  on  either  side  of  the  bend  have  identical 
magnitudes,  thereby  enrorcing  continuity  of  the  current 
flowing  around  the  bend.  Thus  the  current  is  expanded  as 


X 


J(p)  = 


I 

n=l 


J„  pn(p) 


(3.5) 
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Figure  3.2.  Pulse  basis  functions  and  natch  points  for  the 
TM  illuminated  bent  strip. 


V^’  ;*■ «^t^'~<r.i^i 


r'" '  5* 


RWS3 


where 
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0  ,  otherwise 
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and  p  -  ss.  The  coordinate  quantities  in  (3.6)  art 
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>t 
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x  x 


(n-N  )  Ac  ,  N  +1  5  n  -  N+l 

X  X 


(3.7) 


where  is  the  number  of  non-zero  current  pulses  (including 
tne  current  at  the  bend)  along  the  x-directed  portion  of  the 
strip  (see  Figure  3.2).  The  pulse  widths  for  x-directed  and 
c-directed  currents  are 


(3.8a) 


Ac  =  - £ -  , 

N  -  S  +1 
x 


(3-8b) 


respectively.  There  is  a  total  of  N  unknown  currents  repre¬ 
sented  by  pulse  functions.  Substitution  of  (3.5)  into  (3.2) 
permits  the  computation  of  the  vector  potential.  To  compute 
the  scalar  potential,  we  use  a  finite  difference  approxi¬ 
mation  of  the  continuity  equation  to  obtain  a  pulse  repre¬ 
sentation  of  the  charge  density  (Figure  3.3): 
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(3.9) 


where 


ana 


7^  ^  1 


P(P)  = 


di 


I 

n=l 


i  -  J  . 
_n _ n-1 

As 


As  = 


ax 
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L  “c 


,  n  -  N 

’  V 


,  n  '  S  +1 
x 


(3.10) 


WP)  =  < 


1  ,  s  £  s  <  s 

n  n-1 


0  ,  otherwise 


In  writing  (3.9),  we  assume  that 


J  =  J 
0  N+l 


(3.11) 


(3.12) 


Equation  (3.1)  Is  enforced  at  the  centers  of  the  current  sub- 
domains  (Figure  3.2)  on  the  straight  segments  of  the  strip. 
Replacing  the  term  involving  the  derivative  of  the  scalar 
potential  by  a  finite  difference  approximation,  we  obtain. 


me 


-  E. 


(xs  )Ax  -  - juA  (xs  )Ax  --  Q  x  (s  + 


x  m 


x  a 


a  =  1,  2,  ...,  Ns“l  (3.13a) 
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_mc  x  * 

■c.  (cs  )Ac  = 
c  m 


-iuA  (cs  )Ac  -  $ 
c  m 


^X)] 

#»x)]  • 


m  =  X  +1,  X  +2,  ...,  X  (3.13b) 
x  x 


At  the  bend,  we  enploy  a  backwards  difference  approximation 
to  the  derivative  of  the  scalar  potential  at  s  =  0  (on  the 
left  side  of  the  bend  in  Figure  3.2)  and  a  corresponding 
forward  difference  approximation  at  s  =  0+  (on  the  right 
side  of  the  bend  in  Figure  3.2): 


-  Ax 


-E^nC(0)  =  -  jU)A^_  (0)  - 


(0)  -  H-x  x 

Ax/2 


,  ( 3 . i 4a) 


inr  -  _ 

•Ec  (0)  =  ~ju)Ac(0)  - 


(«¥)- 


$(0) 


Ac  /2 


(3.14b) 


Only  one  equation  is  required  at  the  bend  so  that  some 
appropriate  combination  of  (3.14a)  and  (3.14b)  must  be  taken 
Heretofore,  we  have  successfully  avoided  evaluating  any 
potentials  due  to  singular  quantities  directly  at  the  point 
or  the  singularity  (see  Section  II  and  reference  [  5]). 
aquations  (3.14a)  and  (3.14b),  however,  seem  to  require  the 
evaluation  of  the  scalar  potential  directly  at  the  bend 
“here  the  charge  density  is,  in  general,  infinite.  To  avoid 
any  dirficulty  with  this  questionable  procedure,  we  merely 
multiply  (3.14a;  by  Ax/2  and  (3.14b)  by  Ac/2  and  add  the 
resulting  equations,  obtaining 


(Einc(0)Ax  +  ElnC(0)ic\  /A  (0) Ax  +  A  (0)Ac 

X  C  5=  -jw1  *  C 


-  $(c  +  0^-x  .  (3.15) 


In  this  equation  for  the  natch  point  at  the  bend,  the 
scalar  potential  at  the  bend  is  never  needed.  The  use  of 
(3.15)  at  the  bend  is  an  important  development  because  it 
ensures  that  (1)  no  potential  functions  are  evaluated  at 
singularities  of  the  corresponding  source  quantities, 

(2)  it  quite  properly  gives  more  weight  to  the  longer  sub- 
domain  on  either  side  of  the  bend,  and  (3)  it  reduces  to 
(3.13)  when  the  bend  angle  is  zero  and  Ax  =  Ac.  It  will 
also  be  demonstrated  in  a  following  subsection  that  the  use 
of  (3.15)  at  the  bend  yields  numerical  results  which  are 
relatively  independent  of  the  subdomain  size. 

The  vector  and  scalar  potential  terms  can  be  written 
more  explicitly  by  introducing  the  vectors 


p  =  xx 

X 


(3.16a) 


Pc  =  cc 


(3.16b) 


and  a  rather  unconventional  but  convenient  definition  for 
the  Kronecker  delta, 

1  ,  if  s  =  t 


s  t 


0  ,  if  s  ^  t 


where  s  and  t  are  variable  subscripts  rather  than  integers. 
The  potentials  are  then  written  as 


X  -1 


p  =  x  or  c  (3.18) 


The  functions  v  and  are  defined  in  (2.9)  and  (2.10). 

px  pc 

Equations  (3.13)  and  (3.15)  together  with  (3.17)  and 
(3.18)  constitute  a  linear  system  of  equations  which  can  be 
succinctly  written  in  matrix  form  as 

Z1  =  V  '  (3.19) 


31 


I  = 


(3.21) 


and  the  eleaents  of  Z  can  be  determined  from  (3.13),  (3.15), 
(3.17),  and  (3  IS).  Solution  of  the  simultaneous  equations 
(3.19)  determines  the  coefficients  in  the  basis  representation 
for  the  current.  Equation  (3.5). 


Effect  of  Placement  of  Edge  Subdomains  on  Xumerical  Solutions 
The  proper  placement  of  the  end  subdomains  where  one 
requires  the  current  to  vanish  is  critical  for  obtaining 
good  convergence  in  numerical  solutions.  It  has  been  found 
repeatedly  for  such  problems,  (i.e.,  scattering  by  wires 


sfca^c  ■  "  '^~<-X-'?-<‘-+-*?~  -.T'C~^-r''i-''-L,-' -  X^-^.-lV^’'':;y^':'*l'i  —  '.£ 

-  —- . , _ „  ....  ■^..■■w*  iw."  ■K.vjqwj- 


and  TE  strips)  that  requiring  the  current  to  vanish  in  the 
region  within  a  half  subdomain  width  (compared  to  the  adja¬ 
cent  subdomain  widths)  of  edges  correctly  models  the  actual 
electrical  length  of  the  structure.  As  an  example.  Figure  3.4 
compares  the  convergence  of  the  center  current  for  a  normally 
illuminated  TE  strip  where  the  current  is  assumed  to  be 
zerj  within  a  half  subdomain  and  within  a  full  subdomain  of 
the  edge.  The  slow  convergence  that  results  when  the  current 
vanishes  within  a  full  subdomain  of  the  edge  is  attributed 
to  the  fact  that  increasing  the  number  of  subdomains  simul¬ 
taneously  appears  to  increase  the  electrical  length  of  the 
strip.  Hence  it  is  the  length,  and  not  the  numerical  solution 
procedure,  which  is  converging  slowly. 
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entire  strip,  while  at  the  higher  frequencies,  the  edge 
behavior  is  confined  more  to  the  region  near  the  edges. 


Effect  of  the  Bend  on  Numerical  Solutions 

The  bend  in  the  strip  causes  the  current  which  flows 
around  the  bend  to  have  an  infinite  slope,  in  general,  at 
the  bend.  In  Appendix  B,  it  is  shovTi  that  the  surface  cur¬ 


rent  in  a 


numerical  solution  is  the  vector  sum  of  the  actual 


L 


% 


t 

c 


surface  currents  which  flow  on  opposite  sides  of  the  strip. 

In  the  neighborhood  of  the  bend,  the  variation  of  current 
near  the  bend  must  approach  that  of  the  current  near  the 
edge  of  an  infinite  w'edge  of  the  same  angle.  Since  the  vari¬ 
ation  of  the  current  is  more  rapid  on  an  exterior  than  on  an 
interior  wedge,  the  exterior  wedge  current  will  tend  to  pro¬ 
vide  the  dominant  behavior  of  the  current  at  the  bend, 
particularly  when  the  interior  is  in  the  shadow  region  of 
the  incident  field.  Thus  the  exterior  current  near  the  bend, 
from  Appendix  A,  varies  as 


cij 


J(s)' 


s-*-0 


+  B  sgn (s) ■ s 


!  Ti+  |  ( 


(3.22) 


where 


sgn (s)  =  < 


-1 


s  >  0 


s  <  0 


Figure  3.5  shows  this  type  of  behavior  for  current  flowing 
around  a  90°  bend.  As  with  the  edge  behavior,  this  condition 
becomes  confined  to  a  more  limited  region  about  the  bend  as 
the  frequency  increases. 

Figure  3.6  shows  the  current  c..  a  90°  bend  for  the  case 
when  the  illumination  is  incident  along  the  symmetry  plane 
between  the  two  symmetrical  halves  of  the  strip.  Because  of 
the  symmetry  in  both  the  geometry  and  the  incident  field, 
the  current  must  be  symmetric  about  the  bend  and  hence  the 


second  term  in  (3.22)  must  have  a  zero  coefficient.  One  can 


Figure  3.5.  Induced  current  on  a  O.IX  strip  bent  90  at  s/X  -  0  and 
Illuminated  by  a  TE  plane  wave  normally  incident  to  the 
section  of  the  strip  s/X  <  0. 


behavior  (crosses).  The  data  are  matched  at  the 
nd  at  the  third  numerical  value  away  from  the  bond 
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is  adequately  sampled.  Figure  3.9  shovs  that  changing  sub¬ 
domain  size  at  the  bend  also  has  no  appreciable  effect. 
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IV.  FORMULATION  OF  NUMERICAL  SOLUTION  FOR 
SCATTERING  BY  A  BENT  PLATE 


In  this  section  integral  equations  are  derived  for  the 
surface  currents  induced  on  a  bent  rectangular  plate  subject 
to  plane  wave  illumination.  The  solution  is  complete  in  the 
sense  that  all  components  of  current  which  can  exist  on  the 
structure  are  properly  accounted  for.  A  numerical  orocedure 
which  incorporates  the  numerical  principles  developed  in 
Sections  II  and  III  is  described . 


Formulation  of  the  Integral  Equations 

In  developing  the  equations  for  the  bent  plate  structure, 
we  refer  to  the  coordinate  system  of  Figure  4.1,  where  s  is 
a  unit  vector  which  is  perpendicular  to  the  bend  and  which 
lies  on  the  plate.  The  flat  rectangular  plate  is  treated 
merely  as  a  special  case  where  the  bend  angle  is  zero. 

The  currents  induced  on  the  scatterer  are  found  via  the 
numerical  solution  of  a  set  of  coupled  integro— differential 
equations.  Referring  to  the  coordinate  svst em  of  Figure  4.1, 
the  coupled  equations  may  be  written  as 


zinc,-,  ~  I 

( r  )  *  s  =  |-J 


juA(r)  -  VO(r)  «s 


(4.1a) 


(r)  =  -juA  (r)  -  —  CCr) 

>  y  <jy  ’ 


(4.1b) 


r  on  the  plate 
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where 


The  quantities  J  and  o  are  the  surface  current  and  the  surface 

charge  densities  on  the  plate,  respectively,  and  are  related 

through  the  equation  of  continuity. 

To  effect  a  numerical  solution,  the  plate  is  divided  into 

a  seoarate  set  of  patches  for  each  current  conoonent,  J  and  <1  , 

s  y 

as  shown  in  Figures  A. 2  and  A. 3.  Note  that  the  current  patches 
for  and  are  staggered  relative  to  each  other.  The  stag¬ 
gering  of  the  subdomains  has  several  advantages.  Full  sub¬ 
domains  are  used  to  represent  currents  parallel  to  an  edge, 
as  in  the  corresponding  TM  strip  problems-  Half  subdomains 
are  then  used  to  represent  the  zero  currents  of  components 
normal  to  an  edge,  as  in  the  corresponding  TE  strip  problems. 
Furthermore,  subdomains  for  current  components  normal  to  the 
bend  are  wrapped  around  the  bend, ensuring  continuity  of  cur¬ 
rent  flowing  around  the  bend,  while  components  parallel  to  the 
bend  are  represented  by  separate  subdomains  on  either  side  of 
the  bend,  ensuring  the  independence  of  currents  on  opposite 
sides  of  the  bend.  The  staggering  of  the  orthogonal  current 


A5 


SuhdomnJn  scheme  for  H-dirccted  currents. 


Subdomain  scheme  for  y-directed  currents. 


subdomains  is  also  important  because  it  minimizes  the  error 
in  the  central  difference  formulas  that  are  used  to  approxi¬ 
mate  the  derivatives  in  (4.1)  and  permits  one  to  use  the  same 
numerical  algorithm  for  edge  subdomains  as  for  interior  sub- 
domains.  Furthermore,  the  method  seems  readily  generalizable 
to  treating  surfaces  of  arbitrary  shape.  Within  each  patch, 
the  current  is  assumed  constant  so  'hat  the  expansions  for  the 
current  components,  referring  to  Figures  4.2  and  4.3,  are 


J  (r)  ~  s 
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MSUNKY  MSUNKS 

‘2  2 

LJS  =  1  KJS  =  1 


tKJS,LJS  s  \  /  \ 
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(4.5a) 


0  ,  otherwise 
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where 


Ay  = 


N-  1 


(4.7a) 


As  = 


MPB-1 


(4.7b) 


As  = 


MFB-1 


(4.7c) 


M  =  MFB  +  MPB-1  . 

The  quantities  S  ,  S  ,  and  Y,  are  dimensions  of  the  plate  along 
x  c  d 

the  x,  c,  and  y  axes,  respectively,  and  M  and  N  are  defined 

in  Figures  4.2  and  4.3.  The  unknown  pulse  current  coefficients 

jKJS ,LJS  jKJ\,LJY  -n  (44)  may  be  taken  to  represent  the 

s  y 

current  at  the  center  of  their  corresponding  subdomains  or 
directly  on  the  bend  and  centered  in  the  y-direction  for  sub- 
domains  which  straddle  the  bend.  The  substitution  of  (4.4) 
into  (4.2)  permits  the  computation  of  the  vector  potential. 

To  compute  the  scalar  potential,  the  surface  charge  density 
is  needed.  It  is  conveniently  computed,  as  in  the  TE  strip 
case,  by  a  finite  difference  approximation  to  the  equation 
of  continuity.  The  result  is 
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(4.8) 


,re  the  components  of  current  normal  to  the  edge  satisiy 


0  ,  LJ  _  M- 1 ,  LJ  _  jKJ  ,  0  _  j.KJ  ,N-1  _  Q 

Js  s  y  y 


(4.9) 


The  pulse  basis  functions  for  the  charge  are  defined  as 
(Figure  4.4) 


Li5> 

,y  -  <  s  5  sy  +  -5*.  KJ  <  MPB-1 

KJ  2  KJ  2 


As  As 

1  <->’ _ ^  <  <  y  +  _c  Mps  <  KJ 

1  •  SKJ  2  '  S  KJ  2 


PKJ,LJ 
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(4.10) 


0  ,  otherwise 

Substitution  of  (4.8)  into  (4.3)  permits  the  scalar  potential 
to  be  computed  in  terms  of  the  currents.  Equations  (4.1a) 


and  (4.1b)  are  enforced  at  the  match  points 


valid  on  the  surface  in  the  c-y  plane,  and 
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valid  everywhere  on  the  plate.  At  the  bend,  we  require  the 
tangential  components  of  electric  field  along  both  the  x- 
and  c-fiirections  to  vanish.  For  these  field  components, 
we  write  the  directional  derivative  of  the  scalar  potential 
along  the  x-direction  as  a  backward  difference  and  that 
along  the  c-direction  as  a  forward  difference  operator: 
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Note  the  appearance  of  the  scalar  potential  term  4>(s_„,  y  ) 
which  is  evaluated  at  the  bend  where  the  cnarge  is  normally 
singular.  As  discussed  in  Section  II,  this  would  normally  be 
expected  to  introduce  numerical  anomalies  when  pulse  functions 
are  used.  However,  we  need  only  one  equation  at  the  bend  and 
we  obtain  it  by  multiplying  the  two  equations  above  by  As^/2 
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and  as^/2,  respectively,  and  adding  the  resulting  equations 


together.  The  term  involving  the  potential  at  the  edge  then 


disappears  altogether: 
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JS  =  1,  MSUXKY  . 


Obviously,  if  Asx=Asc>  (4.12b)  and  (4.i2d)  reduce  to  (4.12a) 


when  the  bend  angle  is  set  to  a  =  0,  reducing  the  structure 


to  a  flat  plate.  In  order  to  define  the  quantities  appearing 


in  equations  (4.12)  more  easily,  it  is  convenient  to  define 


two  new  vectors 


r  =  xx  -r  yy 


(4.13a) 


r  =  cc  +  vv  , 
c 


(4.13b) 


so  that  r  defines  a  point  on  the  portion  of  the  plate  in  the 


x-y  plane,  and  r^  defines  a  point  on  the  portion  of  the  plate 


in  the  c-v  plane.  Then  the  terns  appearing  in  equations 


(4.12)  are  defined  in  (4.14)  and  (4.1 5)  which  follow. 
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of  the  Kronecker  delta. 
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and  the  *j  functions  are  defined  as 
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p  -  x  or  c. 
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All  the  quantities  appearing  on  the  right-hand  sides  of 
equations  (4.12)  are  now  defined.  The  incident  field  is  ex¬ 


pressed  as 


mine  2  ,  _  ?.  j kn • r 

r.  =  (t0o  +  E^v)  e 


(4.19) 


where  ?  and  <•>  are  the  standard  spherical  coordinates  with 
respect  to  the  x.  v,  z  coordinate  system.  The  components  of 


che  incident  field  in  the  x-y  plane  can  be  found  by  a  simple 
coordinate  transformation  to  be 
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while  the  corresponding  fields  in  the  c-y  plane  are 
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(4.20b) 

Equations  (4.12)  constitute  a  system  of  linear  equations  for 
the  unknown  current  coefficients  in  (4.4). 
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The  system  can  be  written  as  a  matrix  equation 

ZI  =  V  (A. 21) 

where  Z  is  the  impedance  matrix,  I  is  a  column  vector  of  the 
unknown  current  coefficients,  and  V  is  a  column  vector  con¬ 
taining  the  dependence  on  the  incident  field.  The  impedance 
matrix  elements  consist  of  appropriate  linear  combinations 
of  integrals  of  the  form  (4.17).  The  evaluation  of  (4.17) 
is  by  two-dimensional  Gaussian  quadrature.  In  the  case  where 

r  is  interior  to  the  patch  over  which  the  integration  is  to 
P 

be  performed,  the  integrand  has  an  integrable  singularity. 

For  a  discussion  of  the  numerical  evaluation  of  (4.17)  in 
this  case,  refer  to  Appendix  C. 

The  staggered  subdomain  scheme  has  the  advantage  of 
avoiding  the  evaluation  of  any  field  components  or  potential 
quantities  at  edges  or  corners  which  would  be  adversely 
affected  by  the  absence  of  the  correct  singularity  in  the 
current  or  charge  expansion.  The  actual  behavior  of  current 
at  an  edge  is  discussed  in  Appendix  A  while  the  behavior  of 
current  near  a  corner  is  discussed  in  Appendix  D.  In 
Appendix  E  a  possible  correction  factor  for  edge  currents 
is  suggested.  More  extensive  convergence  studies,  however, 
are  necessary  to  validate  the  legitimacy  of  the  formulas 
given  there. 

An  alternative  form  of  Equations  (4.12)  is  given  in 
Appendix  F.  This  form  describes  the  scattered  field  entirely 
in  terms  of  vector  potential  quantities  and  has  advantages 
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NUMERICAL  RESULTS  FOR  PLATE  STRUCTURES 


v 

In  this  section  we  present  numerical  results  for  flat 
and  bent  structures  illuminated  by  plane  waves.  The  numer¬ 
ical  formulation  and  approach  of  Section  IV  was  used  in  the 
computation  of  all  results  shown  in  this  section.  Other 
numerical  approaches  were  tried  during  this  study  but,  for 
brevity,  none  of  those  results  are  given  here  even  though 
they  are  quite  helpful  in  illustrating  the  advantages  of  the 
adopted  approach. 

In  plate  problems,  of  course,  one  is  severely  constrained 
by  the  number  of  subdomains  that  can  be  used  because  of  com¬ 
puter  central  core  limitations.  Thus,  even  though  the  prob¬ 
lem  may  be  numerically  formulated  correctly  in  the  sense 
that  it  will  converge  to  the  correct  answer  as  the  number  of 
subdomains  tends  to  infinity,  in  a  practical  sense,  it  is  a 
high  rate  of  convergence  that  is  important  so  that  one  may 
obtain  "converged"  answers  with  as  few  subdomains  as  possible. 
Figure  5.1  shews  the  convergence  of  the  dominant  component 
of  the  current  at  the  center  of  a  square  plate  at  two  fre¬ 
quencies  and  for  two  angles  of  incidence.  The  good  con¬ 
vergence  is  attributed  primarily  to  the  more  accurate 
modeling  of  the  electrical  size  of  the  plate  that  is  af¬ 
forded  by  the  staggered  subdomain  scheme.  At  the  lower 
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frequency,  convergence  is  achieved  very  quickly;  at  the 
higher  frequency,  the  figure  indicates  that  about  seven 
subdomains  per  wave  1  eng th  are  needed  for  normal  incidence. 

As  the  angle  of  incidence  increases,  more  pulses  are  needed 
net  only  to  represent  the  more  complex  current  distribution, 
but  also  to  adequately  sample  the  incident  field. 

Figures  5.2  -  5.25  show  computed  current  distributions 
for  flat  and  bent  plates  for  two  frequencies  and  for  various 
angles  of  incidence.  In  all  the  figures,  if  the  plate  is 
thought  of  as  folded  out  flat  (the  resulting  rectangular 
plate  outline  is  shown  in  the  figures),  the  value  of  the 
current  at  each  point  is  the  vertical  distance  from  the 
plane  of  the  flattened  plate  to  the  surface  that  is  plotted 
in  the  figure.  The  perspective  view  has  been  chosen  to  ob¬ 
tain  the  best  overall  view  of  each  current  distribution. 

While  use  of  a  digital  plotter  is  mandatory  in  graphically 
displaying  such  a  large  volume  of  data,  this  mechanization 
is  not  without  some  peculiarities  which  require  our  expla¬ 
nation.  First,  the  bend  is  denoted  merely  by  the  notation 
"s  =  0"  which  locates  the  point  at  which  the  bend  (if  present) 
of  the  unfolded  plate  intersects  the  edge  of  the  plate. 

This  intersected  edge  (in  the  left  foreground  of  the  flat¬ 
tened  plate  in  all  the  figures)  is  parallel  tc  the  s-axis 
(Figure  4.1);  the  edge  in  the  right  foreground  in  all  the 
figures  is  parallel  to  the  y-axis.  The  incident  magnetic 
field  intensity  in  all  cases  is  unity  and  each  figure  shows 
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the  largest  positive  and/or  negative  value  of  the  computed 
current  on  the  vertical  scale(s)  on  the  left  hand  side  of 
the  figure.  Note  that  the  plots  of  current  distribution 
stop  half  a  subdomain  from  any  edge  where  they  would  tend 
to  be  singular  if  continued  to  the  edge  and  that  all  com¬ 
puted  values  are  linked  to  computed  currents  at  adjacent 
points  by  straight  lines;  the  intersections  of  the  lines  in 
the  figures  thus  are  points  at  which  the  current  is  computed 
Unfortunately,  currents  which  flow  parallel  to  a  bend  are 
singular  on  both  sides  of  the  bend,  but  the  digital  plotter 
merely  connects  currents  on  either  side  of  the  bend  by 
straight  lines.  If  this  point  is  overlooked,  some  figures 
appear  to  represent  currents  which  are  completely  unphysical 
(e.g.,  see  Figure  11(a)). 

In  comparing  the  dominant  component  of  the  plate  cur¬ 
rent  (i.e.,  the  current  component  parallel  to  the  incident 
electric  field  vector)  to  the  currents  resulting  from  TE 
and  TM  illumination  of  strips  of  the  corresponding  bend 
angle  and  electrical  size,  it  was  found  that  the  plate  cur¬ 
rents  distributions  along  the  plate  center  lines,  i.e., 
along  the  x- ,  c-,  and  y-axes,  are  remarkably  similar  to  those 
of  the  corresponding  TE  and  TM  illuminated  strips.  This 
result  provided  a  convenient  checkpoint  for  assessing  the 
validity  of  many  of  the  computations  of  plate  currents. 

Finally,  we  note  that  some  apparent  anomalies  :'n  the 
figures,  such  as  the  sudden  downturn  of  the  current  normal 
to  the  edges  in  Figure  5.6(b),  have  been  verified  to  be 
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correct  by  comparing  the  current  distribution  to  the 
sponding  TE  strip  problem  (see  Figure  3.8).  When  such 
anomalies  occur  near  corners  (e.g.,  Figures  5.8(a)  and 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report,  we  have  developed  a  numerical  procedure 
for  solving  thin  bent  plate  scattering  problems  using  a 
method  of  moments  solution  of  the  electric  field  integral 
equation.  The  resulting  technique  is  simple,  efficient,  and 
highly  convergent.  Simplicit  is  attained  because  only 
simple  pulse  representations  of  the  current  and  charge  are 
used,  all  derivatives  are  easily  determined  using  finite 
difference  approximations,  continuity  of  current  flowing 
around  a  bend  is  automatically  incorporated,  and  no  special 
condition  on  the  current  is  required  near  edges  or  corners. 
Furthermore,  all  the  integrals  required  involve  only  the 
free  space  Green's  function  in  the  integrand.  Some  effi¬ 
ciency  is  gained  because  of  the  simple  integrals  that  must 
be  performed.  In  addition,  for  flat  plate  structures  or 
for  flat  portions  of  a  bent  plate,  many  of  the  integrals 
that  must  be  computed  may  be  reused  possibly  a  hundred  or 
more  times  in  the  moment  matrix  filling  operation.  Finally, 
some  of  the  efficiency  results  from  the  high  rate  of  con¬ 
vergence  of  the  method  since  this  reduces  the  amount  of 
central  core  and/or  computation  time  needed  to  obtain  a 
converged  solution.  The  authors  expect  that  further  savings 


in  computation  time  should  be  possible  by  replacing  the 


Some  of  the  important  features  of  the  method  which 
should  be  kept  in  mind  in  any  possible  generalizations  of 
the  method  to  arbitrary  structures  are  the  following: 

1-  Currents  normal  to  an  edge  should  be  represented 
by  half  patches.  The  current  is  taken  to  be 
zero  throughout  the  half  patch. 

2.  Currents  parallel  to  an  edge  should  be  repre¬ 
sented  by  a  full  patch.  The  point  inside  the 
subdomain  to  which  the  numerically  determined 
current  is  assumed  to  correspond  should  be 
near  (or  at)  the  patch  center,  away  from  the 
edge,  so  that  the  patch  current  coefficient 
represents  a  finite  current. 

3.  Currents  normal  to  a  bend  should  be  represented 
by  patches  wnich  wrap  around  the  bend,  thus 
ensuring  continuity  of  current  flowing  around 
the  bend.  The  patch  current  coefficient  for 
each  patch  can  be  associated  with  the  value 

of  current  directly  at  the  bend. 

4.  Charge  subdomains  adjacent  to  edges  o.  ’-ends 
are  chosen  such  that  their  boundaries  coincide 
with  the  edge  or  bend.  The  charge  is  computed 
by  a  finite  difference  approximation  to  the 
divergence  of  the  current  entering  and  leaving 
the  charge  subdomain. 

5.  The  moment  matrix  equation  is  derived  by  forcing 

—  s  — 

the  component  of  scattered  field  E  =  -jwA-9<> 
at  the  center  of  and  along  the  direction  of 
the  assumed  current  in  each  subdomain  to  equal 
the  negative  of  the  same  component  of  the 
incident  field  there.  The  directional  derivative 
of  the  scalar  potential  there  is  approximated 
by  an  appropriate  finite  difference  of  the 
scalar  potential  computed  at  the  centers  of 
adjacent  charge  subdomains. 

Several  approaches  for  implementing  these  techniques 
into  a  general  surface  modeling  algorithm  have  been  considered. 
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However,  such  an  approach  is  still  in  a  formative  stage  of 
development  at  this  point  and  is  therefore  not  reported 
here-  This  approach  is  promising  in  that  it  may  be  possible 
to  replace  present  wire  grid  surface  modeling  codes  by  a 
surface  patch  model  code  which  uses  the  electric  field 
integral  equation. 

The  symmetry  of  the  structure  has  not  been  utilized 
in  the  computer  code  developed  from  the  algorithm  described 
in  this  report.  For  the  flat  plate  structure,  there  exist 
two  symmetry  planes  (for  the  bent  plate,  only  one  symmetry 
plane  exists)  which  can  be  used  to  reduce  the  size  of  the 
matrix  that  must  be  stored  in  the  computer.  For  arbitrary 
angles  of  incidence,  four  separate  problems  then  need  to 
be  solved  to  obtain  the  total  solution.  Since  each  sub¬ 
problem  requires  the  same  matrix  fill  time  as  the  original 
problem  without  utilizing  symmetry,  one  merely  trades  re¬ 
duced  core  requirements  for  increased  computation  time. 
Another  technique  for  reducing  storage  limitations  is  to 
use  matrix  partitioning  to  reduce  the  size  of  matrix  that 
must  reside  in  core  at  any  one  time.  Unfortunately,  this 
technique  also  increases  computation  time.  For  larger 
plates  or  better  resolution  or  convergence  of  the  currents, 
such  techniques  may  be  necessary,  however.  On  the  other 
hand,  one  should  be  aware  that  the  extra  effort  that  must 
go  into  implementing  these  features,  grins  very  little 
additional  capability  in  handling  larger  structures-  This 
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fact  is  due  to  the  so-called  "curse  of  dimensions"  that 
always  plagues  such  multi-dimensional  problems. 

Other  extensions  which  can  be  considered  are  methods 
for  treating  a  wire  attached  to  a  plate  structure  in  the 
central  portion  of  a  plate,  at  an  edge,  or  at  a  corner.  A 
code  is  presently  available  for  treating  a  disconnected 
wire  and  plate  configuration  (Appendix  G)  and  it  remains 
only  to  develop  a  technique  for  connecting  the  two  structures 

There  appears  to  be  no  difficulty  in  extending  the 
present  technique  to  treat  a  corner  formed  by  three  mutually 
perpendicular  conducting  planes.  Besides  the  practical 
interest  in  such  structures,  such  a  geometry  would  be  use¬ 
ful  in  numerically  determining  a  "corner  condition",  since 
no  separable  coordinate  system  is  available  for  analyzing 
such  an  exterior  corner. 

In  conclusion,  the  data  given  here  for  current  dis¬ 
tributions  on  a  bent  plate  is  extremely  useful  since  so 
little  information  is  currently  available  as  to  what  happens 
to  current  flowing  on  a  bent  rectangular  plate.  However, 
possibly  more  important  for  future  numerical  wcrk  is  the 
ability  of  the  staggered  current  subdomain  scheme  to  re¬ 
solve  a  number  of  apparently  conflicting  req i.irements  on 
the  current  representation,  thus  leading  to  a  simpler  and 
more  efficient  numerical  scheme. 
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APPENDIX  A 


BEHAVIOR  OF  CURRENTS  NEAR  BENDS  AND  EDGES 


1 


In  this  section,  the  behavior  of  currents  flowing  near 
edges  and  bends  on  conducting  surfaces  is  summarized.  The 
correct  behavior  of  the  currents  is  rigorously  obtained  by- 
requiring  that  the  total  energy  be  bounded  in  any  neighbor¬ 
hood  of  the  edge  not  containing  sources.  We  follow  here  the 
argument  of  Jones,  which,  in  his  words,  is  "plausible  though 
lacking  the  rigor  of  ...  (other)  ...  investigations"  [A  ]. 
Jones  argues  that  the  behavior  of  the  current  near  an  edge 
or  bend  car  be  determined  from  an  infinite  cylindrical  wedge 
of  the  same  wedge  angle  as  the  bend  or  edge.  Figure  A.l 
shows  the  cross-section  of  such  an  infinite  cylindrical  wedge. 
In  the  region  0  i  6  5  2»-3,  the  fields  can  be  decomposed  into 
parts  either  TE  or  TM  to  z  and  can  be  expressed  in  terms  of 


corresponding  potential  functions  which  are  source  free  solu¬ 
tions  of  the  wave  equation  in  cylindrical  coordinates  [  6): 
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The  Bessel  functions  of  the  first  kind  are  used  in  (A.l)  and 
(A. 2)  because  only  they  satisfy  the  finite  energy  condition 
in  the  neighborhood  of  p  =  0  [  *].  The  field  components  are 
then  given  by 
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The  separation  constant  v  and  the  4>-dependent  harmonic 

n 

functions  in  (A.l)  and  (A. 2)  are  chosen  such  that  (A. 5a)  and 
(A.  5c)  both  vanish  at  <■>  =  0  and  4>  =  2tt-S  as  required. 

Referring  to  Figure  A. 2,  we  note  that  a  bend  in  a  plate 
(strip)  forms  both  an  interior  and  an  exterior  wedge.  There 
are  two  components  of  current,  one  flowing  parallel  and  one 
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flowing  perpendicular  to  the  bend.  If  one  notes  that  the 


surface  current  density  is  defined  as  J  =  n  *  H  where  n  is 
a  urit  outward  normal  to  the  plate  (strip),  and  that  for 
small  values  of  p. 
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(A. 6) 


then  one  easily  deduces  from  (A. 5),  using  the  conventions  cf 
Figure  A. 2,  the  following  expansions  for  the  current  near 
the  wedge  tip  at  s  =  0: 
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The  coefficients  in  (A. 7),  of  course,  are  functions  of  po¬ 
sition  along  the  bend.  One  should  keep  in  mind  that  in  a 
numerical  solution,  the  unknown  current  that  is  numerically 
determined  in  a  plate  or  strip  problem  is  actually  the  sum 
of  the  interior  and  exterior  currents.  Except  for  the  con¬ 
stant  terras,  the  exterior  current  typically  dominates  the 
interior  current  for  small  s  because  it  involves  lower  powers 
of  s  or,  in  the  case  of  J^,  is  singular. 

When  a  =  77 ,  the  tip  of  the  wedge  becomes  an  edge  on  a 
vanishingly  thin  slab  of  conducting  material  (Figure  A. 3). 


In  computing  scattered  fields  from  the  currents  flowing  on 
a  metallic  scaiterer,  we  use  the  equivalence  principle  [  6  ] 
and  remove  the  metal,  leaving  the  currents  to  radiate  in 
free  space.  For  the  vanishingly  thin  scatterer,  the  currents 
on  either  side  of  the  edge  then  radiate  from  the  same  surface 
in  space  and  hence  may  be  replaced  by  a  net  equivalent  cur¬ 
rent  which  is  the  appropriate  vector  sum  of  the  currents 
on  either  side  of  the  edge.  Referring  to  Figure  A-4  for 
the  assumed  current  directions  and  using  (A. 7)  with  a  =  r , 
we  obtain 
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=  -2Bis^  -  2B3S  -  ...  ,  s  >  0  .  (A. 9b) 

Thus  the  component  of  (net)  current  parallel  to  the  edge  has 
a  square  root  singularity  while  the  component  of  (net)  cur¬ 
rent  perpendicular  to  the  edge  vanishes  as  the  square  root 
of  distance  from  the  edge. 
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APPENDIX  B 

SURFACE  CURRENT  ON  A  THIN  SCATTERER 


Consider  currents  which  exist  on  a  thin  (maximum  thick¬ 
ness,  At)  scatterer  as  shown  in  Figure  B.l(a).  The  scatterer 
is  assumed  to  have  two  readily  identifiable  "sides”.  Quan¬ 
tities  associated  with  one  of  the  sides  are  denoted  by  a 
positive  sign  superscript,  the  other  side  by  a  negative  sign 
superscript.  Thus,  for  example,  the  surface  currents  are 
given  hy 

—  +  r\  t-  _  + 

J  =  n~  x  H'  <,B.l) 


where 


H~  =  H‘ 


-me 
+  H 


(B-2) 


■“  i  Tl  c 

The  incident  magnetic  field  intensity,  li  ,  is  essentially 

—  s  - 

the  same  on  both  sides  of  the  strip  if  Ax  is  small.  H 
is  the  scattered  field  on  the  positive  (negative)  side  of 

„  -f 

the  scatterer  and  n_  is  a  unit  vector  normal  to  the  scatterer. 
According  to  the  equivalence  principle  [  6  ),  as  far  as  the 
region  exterior  to  the  scatterer  is  concerned,  the  scatterer 
may  be  removed  and  the  surface  currents  J”  placed  on  the 
surface  vacated  by  the  scatterer  (Figure  3.1(b)).  With  the 
scatterer  removed,  these  currents  radiate  in  free  space  and, 
together  with  the  incident  fields,  produce  the  correct  fields 
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(d) 


Figure  3.1.  Sequence  of  steps  to  arrive  at  a  net  equivalent 

current:  (a)  original  scatterer  of  maximum  thickness 

£t,  (o)  scatterer  replaced  by  equivalent  currents 
radiating  in  free  space,  (c)  equivalent  currents  for 
a  vanishingly  thin  scatterer,  and  (d)  replacement  of 
currents  on  either  side  by  a  net  equivalent  current 
representing  the  scatterer. 


as  i 


everywhere  exterior  to  the  scatterer  (and  zero  fields 
interior  to  the  scatterer,  but  this  is  of  no  consequence  in 
the  fo 1 JOwing) . 

Now  in  many  problems,  the  maximum  thickness.  At,  of 
the  scatterer  is  so  small  that  it  is  convenient  to  assume 
that  At  =*  0.  Yhe  result  is  that  the  equivalent  surface 
currents  J  and  J  now  exist  on  the  same  surface  (Figure 
B.l(c))  cannot  be  separately  solved  for  in  a  numerical 

solution  procedure,  (Note  that  when  we  require  the  tangential 
electric  field  to  vanish  at  a  point  on  the  scatterer,  there 
is  no  longer  any  mechanism  to  distinguish  which  side  of  the 
scattetet  we  are  on.)  Instead,  one  solves  for  the  net 
equivalent  current  (Figure  B.l(d)) 

3  =  J+  +  J-  .  (5.3) 

If  one  subsequently  desires  to  know  the  current  on  one  side 
or  the  c^her,  he  must  first  compute  the  scattered  field  from 

HS  =  -  V  *  A  (B.  4) 

l! 

where 


A  = 


J/3<) 


-jk|r-r’ i 


-  dS’ 


(B.3) 


r-r 


scatterer 


_+ 

and  then  form  J  or  J  using  (B.l)  and  (B.2).  Some  care 
must  he  taken  in  the  evaluation  of  Hs  on  the  scatterer  sur¬ 
face;  cn^  should  see,  for  example,  [  7  ]  in  this  regard. 


The  situation  is  somewhat  simpler  if  the  scatterer  is 


planar,  such  as  a  flat  strip  or  flat  plate;  by  symmetry,  one 
has  in  this  cas e 

n*  x  hSZ  =  |  ,  (B . 6 ) 

and  hence  by  ^B.l)  and  (B.2), 


_ -r 

J" 


n 


-inc 


x  H 


(B .  7) 


The  differences  in  the  interpretation  of  the  equivalent 
currents  represented  by  (B.l)  and  (B.3)  should  be  particular¬ 
ly  kept  in  mind  in  comparing  experimentally  measured  currents 
with  numerically  computed  currents.  The  latter  is  of  the 
type  (B.3)  whereas  the  former  is  generally  of  the  type  (B.l). 


APPENDIX  C 


INTEGRATION  OF  THE  SELF  PATCH 

The  integrands  of  the  functions  v  and  ^  are  sin- 

xx  cc 

gular  whenever  the  field  point  is  located  within  the  source 
patch  region.  Bv  a  simple  translation  of  the  origin  these 
self  term  integrals  can  always  be  written  in  the  following 
form : 


(C.l) 


where  As  =  Ax  or  As  =  Ac  depending  on  whether  (C.l)  is  the 


self  term  for  i!>  or  o  .  respectively.  Since  the  integrand 
xx  c  c 


is  symmetric,  y 

self 


mav  also  be  written  as 


ill 


self 


M  4Z  f  7  3  *S 

J2  J2  e-jki-s“  +  yJ 


r? 

k~  +  y  J 


dyds 


(C.  2) 


0  0 

Subtraction  and  addition  of  tht  term 


vields 


self 


=  4 


o  o  L- 


-jk 


[s2  +  y2J 


-  1 


[s2  +  y2]  [f2  +  y2] 


4dyds  • 

(C.  3) 
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The  first  term  in  the  integrand  of  (C.3)  is  now  bounded 
and  smooth,  and  can  therefore  be  integrated  accurately  by 
numerical  quadrature.  The  second  term  in  the  integrand  is 
singular  at  s  =  y  =  0.  Thus  we  can  write  *  Jf  as  the  sum 


of  two  integrals 


>J)  =  li>k  +  i/® 

self  se 1 f  sell 


(C.4) 


where  t*;b  has  a  bounded  integrand  and  .  ,  has  a  singular 
sel f  seii 

integrand . 


We  consider  now  4*^ e ^ f  which  is  given  by 


As  Ay 

r2  r  2 


self 


Urtr 


dyds 


(C.5) 


Integration  with  respect  to  y  yields 


3^ — —  -r  — 


3.'C? 


_  _  _  „  _  __  _ _ _  II I  „  .  « !i.i,u',<  M«jji'iWtBa!»MSBiB8@9* 


Integration  with  respect  to  s  now  yields 
V  ..  =  |4s  1( 


or 


self 


-+-  4 


og{^+[s2  +  (^^. 
(f)  108  |s  +  [=2  +  (f)  ] 

-,S  =  M 

Js  =  0 


-  4s  -  4s  log  s  +  4s 


(C.  8) 


’P  =  2As  log 


self 


+  2  Ay  1 


[f + M +  (€]} 

°8  + [W2  +  W] } 

-  2As  lo,(^)  -  2Ay  log^^ 

•*  {(s)  •  [■  *  («)] } 

••  {(S)  *  [■  *  («)tl  • 


=  2AS  1 


+  2 Ay  log  <  — 


(C .  9} 


Thus,  upon  obtaining  a  numerically  determined  value  for 

,,,  the  value  of  i  mav  be  determined  from  (C.4)  and 

self  selt 

(C. 9) - 
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APPENDIX  l' 
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BEHAVIOR  OF  SURFACE  CURRENT  NEAR  A  CORNER 


In  a  numerica'  solution  to  a  plate  scattering  problem, 
one  might  be  properly  concerned  with  whether  or  not  the 
numerical  solution  represents  (or  is  '  ven  capable  of  repre¬ 
senting)  the  correct  behavior  of  the  current  near  a  corner 
where  two  edges  meet  perpendicularly.  The  corner  current 
has  been  examined  durl.-g  this  studv  and  in  this  section,  we 
show  that  although  no  explicit  condition  on  the  current  is 
imposed  at  the  corners,  the  numerical  solution  does  seem  to 
exhibit  the  correct  behavior  at  the  corners.  The  behavior 
of  the  current  at  the  corners  or  the  "corner  condition"  may¬ 
be  determined  from  the  solution  of  the  problem  of  scattering 
by  a  quarter  plane  {  8  ]  formed  by  taking  one  quadrant  of  an 
infinite  ground  screen.  This  problem  may  be  solved  by  sep¬ 
aration  of  variables  in  spnero-conal  coordinates  [  8  ]  and 
the  eigenfunctions  of  the  resulting  vector  Helmholtz  equation 


can  he  used  to  determine  the  corner  condition  in  the  same 
wav  as  the  cylindrical  wave f unc t ions  are  used  to  determine 


I 


edge  condit ions  on  a  wedge  (Appendix  A). 

Figure  D.  1  shows  a  quarter  plane  lying  between  the 
positive  x-  and  y-axes  in  the  z  =  0  plane.  Points  on  the 
quarter  plane  are  located  by  coordinates  ( C: , C )  where  c  is 


Figure  D.l.  Spnero-conal  coordinates  on  an  infinite  quarter  plane 


Che  distance  to  the  origin  and  6  is  the  angle  that  an 
elliptic  conical  surface  makes  with  the  half  line  y  =  -x, 
z  =  0,  x  <  0  (see  Figure  D.l).  For  points  on  the  quarter 
plane,  the  relationships  between  the  two  coordinate  systems 
are 


(D.la) 
(D. lb) 

(D.lc) 


(The  coordinate  system  of  [  8  }  has  been  rotated  here  by  v/A 
radians  to  better  fit  the  present  discussion.)  According 
to  Satterwhite  and  Kouyoumjian  [  8  ],  the  dominant  behavior 
of  the  surface  current  near  the  corner  is  given  by 


J(r)  = 


-0.186 


6  +  1.228 


i  2 

V 1  4-  sin  6 
sin4> 


*;2W=j 

(D.2) 


where  the  functions  6  „(6)  and  its  derivative  are  periodic 

o  Z 

in  6  and  hence  have  Fourier  series  representations; 


*o2<*> 


($) 


0.96A  sind  +  0.010  sin  34*  +  0.001  sin  5$ 

(D- 3a) 


0.96A  cos6  +  0.030  cos  3$  +  0.005  cos  56 

(D.3b) 
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As  a  check  on  the  consistency  of  (D.2),  we  demonstrate  that 
it  satisfies  the  edge  conditions.  Because  of  the  symmetry 
of  the  expression,  it  is  sufficient  to  examine  the  edge 
along  the  y  axis  where  $  =  0 .  For  small  <•> , 


$->-0 


i>-*0 

2 

2 


Hence , 


*  -0.186,  , 

A  C  |  —  4>(J>  + 


$-♦■0 


[-*  ♦  ^  5] 


=  A  p 


-0.186 


[  -  $  +  1-228  /lf  p]-  (D*A) 


We  see  that  the  6  component  (normal  to  the  edge)  vanishes 
as  the  square  root  of  x  while  the  p  component  (parallel  to 
the  edge)  varies  as  the  reciprocal  of  the  square  root  of  the 
distance  from  the  edge,  in  agreement  with  the  edge  conditions 
(Appendix  A) . 

In  order  to  compare  with  the  numerical  results  obtained 
in  Section  V,  we  examine  the  y-component  of  J.  Noting  that 


o  = 


XX  +  vv 


{ 


(D. 5a) 


2  2 
x  +  v 


6  = 


_  yx  -  xy 


{■ 


2  2 
x  +  y 


(D. 5b) 
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one  readily  determines  from  (D.2)  that 


m 


2  2 

•1^  =  A  (x  +  y  ) 


where 


-0.593 


x  *  . 
L  o. 


(4>)  +  1.228 


(0-6) 


ft 

I 


0  =  cos 


r  2  2 

x  +  y 


In  order  to  get  a  feeling  of  the  current  shape  in  the  neigh¬ 
borhood  of  the  corner,  (D.6)  is  plotted  in  Figure  D.2.  It 
is  interesting  to  note  in  the  figure,  that  indeed  the  current 
vanishes  as  the  square  root  of  distance  from  the  x-axis  and 
has  a  square  root  singularity  on  the  y  axis,  but  also  that 
along  a  radial  line  approaching  the  corner,  the  current  is 
singular  as  o  Comparison  with  some  of  the  low  fre¬ 

quency  plate  results,  say  Figure  5.2(a)  shows  that  there  is 
at  least  qualitative  agreement  between  the  computed  corner 
behavior  and  the  theoretical  corner  condition. 


'•S* 


ehavior  of  the  magnitude  of  the  y-component  of  current  (j 
car  the  corner  of  an  infinite  quarter  plane  located  in  the 
irst  quadrant  of  the  x-y  plana. 


APPENDIX  E 

PATCH  CURRENT  CORRECTION  FACTOR  FOR  CURRENTS 
FLOWING  PARALLEL  TC-  \N  EDGE 


ap 

:i 


In  f  5  ],  it  was  empirically  found  that  for  TK  scattering 
by  infinite  strips,  the  coefficients  of  a  pulse  expansion 
of  the  current  converge  to  the  correct  current  evaluated  at 
the  center  of  the  pulse,  except  for  the  subdomains  near  an 
edge.  In  the  edge  subdomain  ,  convergence  to  the  correct 
current  is  achieved  only  when  the  correct  singular  edge  be¬ 
havior  of  the  current  is  explicitly  incorporated  in  the 
current  representation  in  these  subdomains.  However,  only 
the  edge  subdomain  currents  appear  to  be  affected  by  the  use 
of  a  more  accurate  representation  there.  Hence,  the  edge 
currents  can  h^i  corrected  a  posteriori  by  means  of  a  cor¬ 
rection  factor  which  is  the  ratio  of  moment  method  impedance 
matrix  elements  for  the  self  terms  for  the  edge  subdomains 
with  and  without  the  correct  singular  edge  current  repre¬ 
sentation.  Later  investigations  have  shown  that  these 
empirically  observed  results  have  some  theoretical  basis, 
as  well. 

It  appears  reasonable  to  expect  that  results  analogous 
to  those  above  hold  for  current  components  parallel  to  the 
edges  of  a  rectangular  plate.  With  this  assumption,  for 


115 


■rat; 


edge  patches,  convergence  to  the  correct  current  would  re¬ 
sult  only  if  the  correct  edge  singularity  is  incorporated 


in  the  current  representation.  If  no  other  currents  are 
affected  by  this  special  treatment  of  the  edge  currents, 
then  there  again  exists  a  correction  factor  for  the  currents 
near  the  edge  which  is  determined  by  a  ratio  of  the  impedance 
matrix  self-terms  for  a  patch  with  and  without  the  correct 
singularity  in  the  representation.  Specifically,  we  must 
compute  the  ratio 


where 
from  a 
and  I 

s 

s ingul 
that  I 
at  an 
is  of 
is  sou 
are  su 


pulse  representation  of  the  singular  edge  current 
represents  the  more  accurate  value  obtained  from  the 
ar  pulse  representation  in  parentheses  in  (E.l).  Note 
n  the  assumed  representation  the  current  is  singular 
edge  located  at  x  =  -Ax/2  but  that  the  singular  pulse 
unit  height  at  x  =  0  where  the  value  of  the  current 
ght.  If  the  frequency  is  low  enough  or  the  patch  sizes 
fficiently  small  such  that 


r 


max  kR 


+  Ay2  <<  1 


^»®:'*':«S»>-y;!S^g^5jjtiv^g!p^^^ES>5: 


f^B^^gggggsgsa^ 


then  the  exponential  function  can  be  approximated  by  the 
first  term  in  its  series  expansion  about  R  =  0  so  that  we 
ob  tain 


_E  % 


Ax  Ay 


•ja?  r 2  r 2 _ i _ 

_4  A  1^T^'1pT7' 


dvdx 


Ax  Ay 

/T/'  ' 
•'Ax  A 

2  2 


•vAx  _*4z  Jx2  +  y2' 

->  9  ¥ 


dydx 


(£-2) 


The  denominator  of  (E.2)  can  be  readily  evaluated  as  shown 
in  Appendix  C  and  is  given  by 


Ax  *22. 
2 


/  / 


rAx 

2 


_4  Jx2  +y~ 

?  * 


dydx  =  2Av  in 


+  2Ax  2-n 


Ax 

Ay 


.4i 

Ax 


+  1 


{8 


+  1 


J(E-3) 


The  integral  appearing  in  the  numerator  of  (E.2)  is  some¬ 
what  more  difficult  to  evaluate.  An  exact  solution  can  be 
found  in  terms  of  an  infinite  series,  but  the  following 
form  is  more  practical  to  evaluate: 
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44  •> 


Ax  A  v 


Ax  ^Ax 
2  ”  2 


dydx  = 


/  .Ax  2  ^  2 

XT~  Vx  +’/ 


=  2  y?  Ax  dn 


+  A^Ax  -  2Ax  dr.  — - -  +  2V2  Ax  dn  1  +Jl+  (jA 

V?  +  ij  [  V  W/ 


Ax 

2  di 


+  Y2A>: 


-  tnl  +1  1+  ~ 


-  dx  . 
(E.4) 


The  only  integral  remaining  on  the  right  hand  side  of  (E.4) 
has  a  non-singular  integrand  and  can  be  accurately  integrated 
numerically.  It  is  possible  to  express  the  ratio  completely 
in  terms  of  the  ratios  of  the  subdomain  dimensions  as 


follows : 


•"IS  I 


V?  dn  +  2{2  -  dn  +  ifT  dn  1  +  J  1  +  (7^) 

h1/i7  ^T1  +V14'j21  -  Xu 

.  Vta.i  iRF 

_£  £  - -  - - - - — - -  - ^ -  . 

Ts  £  *- [l +  [t5 +  #"-)  <E-5> 

Usir.g  (E.5),  one  can  determine  the  ratio  1  /I  as  a  function 

p  s 

of  Ax/Ay,  as  shown  in  Figure  E.l.  The  result  obtained  in 
this  manner  is,  of  course,  a  quasi-static  approximation  since 
only  the  first  term  in  the  series  expansion  of  the  free 


_£  $ 


Figure  E.l.  Plot  of  the  quasi-static  patch  current  correction  factor  for 
currents  flowing  parallel  to  an  edge  as  a  function  of  the 


space  Green's  function  is  used-  No  attempt  has  been  made 
to  verify  that  (E-5)  does  indeed  predict  the  correct  edge 
current  ncr  has  the  approach  been  extended  to  correct  cur- 


£ 


i- . 
* 


rents  flowing  parallel  to  a  bend, 
be  worthwhile  to  extend  the  capabil 
described  herein  for  treating  plate 


Such  investigations  would 
ities  of  the  techniques 
structures . 
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APPENDIX  F 

ALTERNATE  FORMS  fOR  THE  INTEGRAL  EQUATION  FOR  A  3ENT  PLATE 


It  can  be  shown  by  means  of  a  series  of  tedious  but 
straightforward  manipulations  that  the  scalar  potential  de¬ 
fined  by  (4.3)  and  (4.8)  jar.  be  written  in  the  alternative 


form 


»(r)  =  -  —  j 

j«ps  ! 


r,xf. .  Asx  xf-  a 

Axlr  "  T~  V  -  \V  -  ~2~  SJ 


Ac(7  +  ~r  a)  ~  Ac(;  - 


,  \(f  *  t  y)  -  %(?  -  ¥  i) 


+  ®X(r)  +  $C(r) 


(F.l) 


X  c 

The  quantities  A  and  A  arise  from  partitioning  the  vector 

x  c  •  ° 

potential  as  follows: 


A  =  AX  +  AC 
XXX 


(F. 2a) 


A  =  AX  +  AC 

c  c  c 


(F. 2b) 


where  the  superscri-ts  denote  partial  potentials  arising 


from  currents  which  flow  only  in  the  direction  indicated  by 
the  superscript.  Thus,  the  vector  potential  A*  (p  =  x  or  c) 
is  defined  by  the  first  two  lines  of  (4.14b)  while  A^  is 
associated  with  the  last  two  lines  of  (4.14b).  The  terms 
$  (r)  and  <>  (r)  are  scalar  potentials  produced  by  currents 
which  flow  around  the  bend  in  the  plate: 


iv 


<*>  "  4^e 


N  1  jMPB- 1 ,  LJ  /•yLJ  +  2 
s 


I 

LJ  =  1 


Li  S 


/LJ  2  f 

«*As 

S  X 


s  iyy 
rLJ  2 


- 1  k  r-  r 


{r-r* 


dx’dy 


(F.3a) 


v_i  vs  Av 

*  A  jMPB-1 ,LJ  p '  L J  2 


t:S 


I 

LJ=  1 


As 


r  ■  l 


.""Ik !  r-r’  j 


dx  ’  dy 


FAs  j  r-r’ | 

s  t\y  _  c  1  1 

'  LJ  ~  2  ~2~  (F.  3b) 


These  terms  are  related  to  the  consequences  of  the  numerical 
approach  on  the  "electrical  completeness"  that  is  always 
implicit  in  the  use  oi  the  Lorentz  gauge  condition  to  replace 
scalar  potential  by  vector  potential  quantities  [  9  ]•  Note 
that  the  moments  of  the  effective  charges  from  which  the 
potentials  (F.3a)  and  (F.3b)  arise  are  equal  in  magnitude 
but  of  opposite  sign  and  hence  cancel  as  the  bend  angle,  a, 
and  the  difference  in  subdomain  lengths,  As^  -  As^,  approach 
zero.  The  contribution  of  these  terms  also  vanishes  when 
As^  and  Asc  both  approach  zero.  Hence  (F.l)  is  a  numerical 
approximation  to  the  Lorentz  gauge  condition 
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**>  (  r  )  =  - 


jOJ^iE 


3A  3Al  3a 


3s 


21  +  __c  +  _v 
3sc  3y 


jwye 


V*  A 


(F.4) 


If  one  writes  (F.l)  as 


0(r)  =  - 


I 


juye 


AA  (  r  ) 
x 


As 


AA‘(r)  AA  (r) 

+  — 7" -  +  — * - 


Ay 


+  0X(r)  +  0C(r) 


(F.5) 


where  the  central  finite  difcerence  operators  in  (F.5)  are 
defined  by  comparing  to  (F.l),  then  Equations  (4.12)  can  be 
written  in  the  alternative  form 
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We  emphasize  that  Equations  (F.6)  are  numerically  equivalent 
to  (4.12).  This  present  form  has  the  advantage  of  illus¬ 
trating  explicitly  the  fact  that  the  scalar  potential  inte¬ 
grals  are  nothing  more  than  vector  potential  integrals 
evaluated  at  different  observation  points.  On  the  other 
hand,  (4.12)  is  a  more  compact  form  and  does  not  require  a 
special  treatment  of  the  edge  currents  as  does  (F.6). 
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APPENDIX  G 


DETERMINATION  OF  CURRENTS  EXCITED  ON  A  WIRE/PLATE  STRUCTURE 


In  this  appendix,  we  investigate  the  problem  of  de¬ 
termining  the  currents  excited  on  a  structure  which  com¬ 
prises  a  straight  wire  and  a  bent  plate  and  which  is 
illuminated  by  a  known  incident  field  (E1 ,  H1) .  Coupled 
integro-di f f erential  equations  are  formulated  for  this 

problem  in  terms  of  the  unknown  surface  current  J  on  the 

s 

plate  and  the  unknown  axial  current  I  on  the  wire-  These 
equations  are  solved  numerically  by  the  methods  devel¬ 
oped  in  this  report  and  sample  results  are  presented. 

The  plate  and  wire  are  depicted  in  Figure  G.l,  which, 
in  addition,  serves  to  define  various  quantities  pertinent 
to  the  present  problem.  The  wire  is  thin  relative  to  its 
length  and  to  the  wavelength  of  the  electromagnetic  field, 
and  the  plate  is  vanishingly  thin.  Both  the  wire  and  the 
plate  are  assumed  to  be  perfect  conductors,  and  the  usual 
approximations  of  thin-wire  analyses  are  employed  here. 

A  Cartesian  coordinate  system  is  located  on  the  plate 
as  shown  in  the  figure  where  one  sees  that  the  bend  is 
coincident  with  the  x  axis.  One  surface  of  the  bent  plate 
in  Figure  G.l  is  labeled  and  is  in  the  x-y  plane  while 
that  labeled  is  declined  by  an  angle  a  from  the  y-axis. 
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Points  (x,y)  on  are  located  in  the  usual  way, 

—  A 

r  =  xx  +  yy,  while  those  on  Sjj  are  located  by 


r  =  xx  +  cc 


(G .  1 ) 


where  c  is  measured  along  the  c  axis  in  the  direction  of 
the  unit  vector  c  defined  bv 


c  =  y  cosa  -  z  sine 


(G .  2  ) 


The  unit  vector  normal  to  the  surface  is  simply  n  =  2 
while  that  normal  to  the  surface  S  is 


n  =  x  *  c 


(G.  3) 


The  geometry  of  the  thin,  cylindrical  wire  is  specified 

by  its  radius  a  and  by  the  location  of  its  lower  and  upper 

endpoints,  (x„,  v  ,  z  )  and  (x,  ,  y,  ,  z,),  respect  ivelv . 
a  3  a  0  D  o 

From  knowledge  of  these  endpoints,  one  can  determine  other 
needed  geometric  quantities: 


LENGTH:  L 


L  ■  [(sb  -  +  (  >'b  -  >a)2  +  (*k  -  za)2]  (G-i3) 


UNIT  VECTOR:  l 


1  s  i  [(*b  -  -.)•  +  (yb  -  >a)> +  (%  -  %)^)  <G-4b) 
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WIRE  CENTER: 


r 

c 


c  '  T  [(Xa  +  *b)X  +  (■’a  +  +  (Za  +  zb)Z]  <G-4c> 


Formulation 

The  total  electric  field  tangential  to  the  surfaces 
of  the  perfectly  conducting  plate  and  wire  must  be  zero. 
This  can  be  expressed  as 

n  *  E  =  0  on  and  (G.5a) 

and 

£*E  =0  on  wire  (G-5b) 


or  as 

n  x  eS  =  -  n  *  E*  on  and  Sjj  (G.6a) 

and 

i*ES  =  -  Z,*E^  on  wire  (G.6b) 

where  E  is  the  total  electric  field,  E1  is  the  specified 

—  s 

incident  field,  and  E  is  the  scattered  electric  field 
produced  by  the  currents  and  charges  induced  cn  the  plate 
and  wire.  The  total  scattered  electric  field  can  be 
written 

ES  =  -juA  -  (G.7) 

where  A  is  the  magnetic  vector  potential  and  $  is  the 
electric  scalar  potential.  A  is  due  to  both  the  wire 
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Attention  is  called  to  the  tact  that  the  plate  is  heat  and, 
thus,  that  Equation  (G.8)  oust  be  enforced  on  both  planar 
surfaces  with  appropriate  interpretation  of  the  unit  nornal 
vector  i  discussed  above.  In  addition,  the  integrals 
over  the  plate  in  Equations  (G.9)  must  be  taken 
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with  the  exact  kernel  K(£-£’>  above  given  by 
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(G.8)  and  (G.ll)  constitute  a  set  of  coupled,  vector 
integro-di f ferential  equations,  which,  in  principle,  can 
be  solved  for  the  unknown  current  density  on  the  plate 
and  axial  current  I  on  th<  wire. 


Satapl  e  Results 

By  a  simple  extension  of  the  numerical  techniques 

developed  elsewhere  in  this  report  for  solving  the  bent- 

plate  equations,  one  can  solve  (G.8)  and  (G.ll)  and  thereby 

determine  J  on  the  plate  and  I  on  the  wire.  Tn  the 
s 

following  paragraphs  are  presented  data  illu;*  -a' ' ve  of 
the  behavior  of  current  on  the  wire  and  plate  xn  selected 
configurations. 

In  Figures  G.2  and  G.3  are  shown,  respec tive ly ,  the 
plate  current  on  a  0.5A  *  0.5*  square  plate  and  a  wire 
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current  is  shown  in  Figure  G. 


of  length  L  =  X  and  radius  a  =  0.001X,  due  to  illumination 
which  is  normally  incident  upon  the  plate  but  which  does 
not  excite  the  thin  wire  directly.  The  wire  is  perpen¬ 
dicular  to  the  plate  and  is  located  above  the  point 
(-0.1  25X  ,  -0.125X)  on  the  plate.  Other  details  of 
plate/wire  geometry  and  incident  illumination  are  recorded 
in  the  labels  of  Figures  G.2  and  G.3.  Since  the  incident 
electric  field  E1  is  normal  to  the  wire,  it  does  not  ex¬ 
cite  the  wire  directly  and,  hence,  the  sole  illumination 
of  the  wire  is  due  to  scattering  from  the  plate.  One 
recalls  that  the  antiresonant  current  on  a  one-wavelength 
thin  wire  is  an  odd  function  and  can  be  induced  by  only  an 
odd-function  excitation  (with  respect  to  the  wire  center)  , 
while  the  forced  response  n  a  wire  of  any  length  is  even 
for  an  even- f unc t ion  excitation  and  is  odd  for  an  odd- 
function  excitation.  Since  the  only  illumination  received 
by  the  wire  is  the  electric  field  scattered  from  the  plate, 
wnich  decays  rapidly  along  the  wire  axis,  one  expects  a 
relatively  significant  odd-function  excitation  compared 
with  the  even-function  part.  In  addition,  the  antiresenant 
(sinusoidal)  current  is  more  readily  excited  on  a  one- 
wavelength  wire  than  is  the  forced  response.  As  ore  would 
anticipate  from  the  brief  discussion  above,  the  current  on 
the  wire  illustrated  in  Figure  G.3  is  essentially  sinusoidal 
The  slight  deviation  from  a  sinusoidal  distribution  is  due 
to  the  even- fun c t ion  portion  of  the  forced  response.  The 
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place  current  depicted  in  Figure  G.2  is  almost  what  it 
would  be  if  the  wire  were  not  present;  small  differences 
can  be  observed  primarily  near  the  wire  location  and  in 
the  small  asymmetry  introduced  by  the  presence  of  the 
wire  . 

With  the  plate/wire  geometry  as  stated  above  but 
with  the  incident  electric  field  now  parallel  to  the  wire 
and  perpendicular  to  the  plate  (see  labels  in  Figures 
G.4  and  G.5),  one  observes  essentially  the  expected 
shifted  cosine  current  on  the  one-wavelength  wire,  which 
is  the  forced  response  due  to  a  constant  illumination. 

The  slight  deviation  from  the  shifted  cosine  manifested 
in  the  asvmmetry  of  the  current  of  Figure  G.5  is  due  to 
the  fact  that  a  small  antiresonant  current  is  induced  by 
the  odd-function,  axial  electric  field  which  exists  on 
the  wire  due  to  secondary  scattering  from  the  plate.  The 
plate  current  (Figure  G-4)  is  excited  only  by  scattering 
of  the  incident  field  from  the  wire  with  no  direct  con¬ 
tribution  from  the  incident  illumination  itself.  Notice 
that  the  plate  current  is  symmetric  about  the  diagonal 
passing  through  (-0.125'  ,  -0.125a)  over  which  the  wire  is 
located  and  that  this  current  is  significantly  smaller, 
as  expected,  than  the  current  depicted  in  Figure  G.2. 

Next  we  consider  the  same  0.5X  *  0.5-1  square  plate 
but  with  a  wire  parallel  to  it.  The  wire  length  and  radius 
are  0.5X  and  0.00i>,  respectively,  its  center  is  0.1X  above 
the  point  (-0.125X  ,  0),  and  its  axis  is  parallel  to  the  v- 
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W.-v 


Induced  on  «  wire  nenr  n  euuare  pinto, 


axis.  The  illumination  is  normally  incident  upon  the 

«  £  A 

plate  with  polarization  such  that  E  *£  =  0;  so  the  only 
wire  excitation  is  due  to  scattering  from  the  plate.  A 
half  waveler-gth  wire  is  resonant  length  and  the  resonant 
current  is  an  even  function  with  respect  to  the  wire  cen¬ 
ter.  This  resonant  current  can  be  excited  only  by  an 
even- func  t  ion  illumination.  Vie  note  from  Figure  G.6  that 
the  presence  of  the  wire  has  little  or  no  discernaDle 
effect  on  the  plate  current.  Also,  the  distribution  of 
plate  current  is  such  that  the  illumination  of  the  wire 
is  an  odd  function  which  can  induce  no  even-function 
resonant  wire  current.  The  wire  current  seen  in  Figure  G.7 
is  entirely  forced  and  is  an  odd  function  about  the  wire 
center  as  expected.  This  non-resonant  current  is  small 
and  should  not  be  viewed  as  sinusoidal  (of  wavelength  X), 
even  chough  its  shape  appears  to  be  (of  wavelength  X/2). 
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The  pinto  current  in  shown  in  Figure  G.6. 
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